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INTRODUCTION 


The YF-12 airplane is a twin-engine airplane with variable-geometry engine 
inlets and Mach 3 cruise capability . Reference 1 reported that difficulties were 
experienced in maintaining airspeed and altitude in high speed flight . The report 
suggested several factors that may have contributed to the piloting difficulties , 
including atmospheric disturbances , airplane/propulsion system interactions , and 
aircraft long-period dynamic characteristics . The effects of an interaction of the 
propulsion system, including the variable-geometry inlets, with the airplane’s 
lateral-directional characteristics are reported in reference 2 . Similarly , the 
effects of the propulsion system's interaction with the airplane’s longitudinal char- 
acteristics needed to be investigated . 

As a result, the NASA Dryden Flight Research Center conducted a series of 
flight tests to investigate the YF-12 airplane’s short-period, phugoid, and height 
modes . (The height mode is an aperiodic mode due to the change in airplane char- 
acteristics with altitude . The change results primarily from the change in atmos- 
pheric density with altitude . ) The tests were conducted at a Mach number of 
approximately 2 .9 and with the propulsion system inlets in the fixed or automatic 
configurations . This report presents the stability and control derivatives that were 
obtained for the velocity and altitude degrees of freedom and the standard short- 
period derivatives . The effect of the inlet’s configuration on the derivatives was 
also determined. The aircraft’s modal characteristics were calculated to examine 
the significance of the derivatives . 


SYMBOLS AND ABBREVIATIONS 


Physical quantities in this report are given in the International System of Units 
(SI) and parenthetically in U .S . Customary Units. The measurements were taken 
in Customary Units . 
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dynamic pressure, N/m (Ib/ft ) 
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reference wing area , m (ft ) 
stability augmentation system 
sampling rate, points/sec 
total observation time , sec 
time, sec 
uncertainty level 
average uncertainty level 
velocity, m/sec (ft/sec) 
airplane gross weight, N (lb) 
weighting function 

^ Force along X-axls (ft/seo^) 
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calculated observation vector 

= Force along Z-axis ^ ^ad/sec^ 
mVo 

= = V’ 

measured observation vector 
angle of attack , deg or rad 
angle of sideslip , deg 
incremental change 

incremental bypass door movement , percent of maximum travel 
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elevator position , deg 


e 

p 

t 

Superscript: 
[ ]* 

Subscript: 

0 


damping ratio 

pitch angle , deg or rad 

3 3 

atmospheric density, kg/m (slug/ ft ) 
height mode time constant , sec 
damped natural frequency, rad/sec 

natural frequency , rad/sec 

matrix transpose 
average condition 


A dot over a quantity denotes the time derivative of that quantity . 


AIRPLANE DESCRIPTION 


The YF-12 aircreift (figs. 1 and 2) is a twin-engine, delta-wing airplane with 
Mach 3 cruise capability. The propulsion system consists of two J58 engines with 
axisymmetric , variable- geometry , mixed-compression inlets. Normal-shock posi- 
tion in the inlet is controlled by the movement of the centerbody spike and the 
bypass doors located on the forward part of the nacelle . As shown in figure 3 , 
spike position and duct pressure ratio are commanded by the inlet computer as 
functions of Mach number , angle of attack and sideslip , and normal acceleration . 
The bypass doors are operated in a closed-loop manner to maintain the commanded 
duct pressure ratio . For short-period maneuvers , such as elevator pulses during 
which Mach number is nearly constant , the bypass doors and spike move together 
primarily as a function of angle of attack . In phugoid maneuvers , where the 
airplane changes primarily in Mach number and altitude , the bypass doors and 
spike move primarily as a function of Mach number . 

Elevens provide pitch and roll control . Elevator deflection was taken as the 
average symmetric deflection of the elevens, and aileron deflection was taken as 
the average differential deflection of the elevens . Yaw control was provided by 
two all-movable vertical tails. A stability augmentation system (SAS) provided 
some artificial stability in pitch and yaw and rate damping in pitch , roll , and yaw . 
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INSTRUMENTATION 


The airplane parameters measured included pitch attitude, pitch angular rate, 
normal and axial linear acceleration at the center of gravity , and elevon and bypass 
door position. The instrumentation was alined with the aircraft body axis (fig. 4) . 
Velocity, altitude, and angle of attack were derived from measurements of pressure 
and temperature at the nose boom . The pressure measurement used to derive angle 
of attack was corrected to the center of gravity location . A time lag correction of 
0 . 4 second was obtained from reference 3 and applied to angle of attack for the 
flight conditions of this study . The repeatability of the parameter measurements 
was within the resolution of the digital recording system . The resolution of this 
system is indicated in the following table: 


Parameter 

Recording system 
resolution 

a, deg 

0.05 

V, m/sec (ft/sec) , . . 

1.0 (3.4) 

q , deg/ sec 

0.04 

9, deg 

0.06 

h, m (ft) 

3.5 (11) 



0.010 

^'3 

0.008 

Sg, deg 

0.025 

^bp ’ 

0.30 


TEST DESCRIPTION 


Three phugoid maneuvers that began at an altitude of approximately 
21,900 meters (72,000 feet) and a Mach number of approximately 2.9 were analyzed. 
The phugoid motion was excited by opening the bypass doors to the full open posi- 
tion for approximately 10 seconds and then returning them to the trim position. 

The drag increment from the bypass door pulse provided phugoid oscillations with 
amplitudes of approximately ±0.02 in Mach number and ±600 meters (±2000 feet) in 
altitude. At first, elevator pulses were used to excite the phugoid motion; however, 
it was difficult to return the elevator to the trim position , and the resulting mistrim 
caused large altitude excursions . The bypass doors were much easier to return to 
the trim position , and using them to excite the phugoid motion resulted in oscilla- 
tions centered about the trim condition . The pitch stability augmentation system 
and inlet control system configurations for the tests are summarized in the follow- 
ing table: 
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Case 


Inlet' control mode 


A 

B 

C 


Pitch SAS 


On 

Off 

Off 


Automatic 
Automatic 
Inlets fixed 


Each phugoid maneuver had between two and three cycles of free oscillation. Roll 
attitude was controlled by the roll autopilot to preclude inadvertent inputs to the 
elevator through lateral stick inputs. For the automatic inlet maneuvers, the bypass 
doors modulated about a scheduled position. For the fixed inlet maneuvers, a 
bypass door position 18 percent above the nominal value was used to provide an 
inlet stability margin during the tests . For all maneuvers , the automatic fu^ 
sequencing system kept the center of gravity approximately constant at 0.24c. 

Eight elevator pulses were performed with the pitch stability augmentation 
system off to determine the airplane’s short-period characteristics. The inlet 
control system configurations for these cases are shown in the following table: 


Case 

Inlet control mode 

Al, A2 

Automatic 

Bl, B2 

Automatic 

Cl, C2 

Inlets fixed 

D, E 

Inlets fixed 


Cases A1 and A2 were performed immediately before phugoid case A and with the 
same inlet control configuration. Similarly, cases B1 and B2 correspond to 
phugoid case B , and cases Cl and C2 correspond to phugoid case C . 


METHOD OF ANALYSIS 


The stability and control derivatives were determined with a maximum likelihood 
estimation method which was a modification of the computer program described in 
reference 4 . The modification consisted of adding the altitude degree of freedom so 
that the altitude derivatives could be determined . Two methods were used to 
evaluate the effects of the inlets . In the first method , the standard stability and 
control derivatives for the fixed inlet configuration were compared to those for the 
automatic inlet configuration. The differences in the derivatives were attributed 
to the inlet . In the second method , the inlet was considered as a separate control 
and control derivatives corresponding to the inlet were obtained from the automatic 
inlet configuration maneuvers . The average bypass door position was used as the 
control parameter representative of the entire inlet system . 
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Because of the large difference in the duration of the short-period and phugoid 
maneuvers , the short-period maneuvers were analyzed separately . A typical short- 
period maneuver took approximately 20 seconds , whereas a phugoid maneuver took 
about 300 seconds. 

The axis system used in the analysis and the positive values of the measure- 
ments are indicated in figure 4. The equations of motion used were as follows: 

“ = V ■ (l^ ®o)® * \p®f>P * ^0 

V = + X^,V - (g oos 0^)6 * X^h + Xg 6^ Xg B t x„ 

e bp ^ 

^ CL V q^ h o e o, bp 0 

e bp ^ 

0 = q + 


h = -VqC + VqG + h 


0 




s + Sx. 
e 8, bp , 
e bp ^ \ 0 




a = — fx a + 
X g a 


= — X a + X-^y + X, h + X„ 6 + X„ 6, 

“ V h 8 e d, bp , 

e bp ^ \ 0 


’] \ 


In addition to the bias terms (Z^, X^, Mq, 0^, h^, , and ), the initial condi- 

tions for a, V, and q were determined. For the short-period maneuvers, the 
derivatives Z , X , M , M , M„ , Z„ , X„ , and M„ were determined, and 
“ “ “ ^bp ®bp % 

the altitude and velocity derivatives were kept fixed at zero. For the phugoid 
maneuvers, the velocity, altitude, and bypass door derivatives were determined 
with the short-period derivatives kept fixed at the average values . 

The estimation method minimized the cost functional 
* 

“ y(0j “ yftjjdt where the observation vectors z and y 
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included a, V , q, Q , h, a , and a . The values of the weighting matrix, D. , are 

TX X X 

shown in the following table for the short-period and phugoid maneuvers: 


Observation 

Short-period maneuvers 

Phugoid maneuvers 

Diagonal values of the matrix 

a, rad 

200,000 

200,000 

y, m/sec (ft/ sec) . . 

0.54 (0.105) 

0.22 (0.02) 

q, rad/ sec 

1,000,000 

2,000,000 

0, rad 

1,000,000 

200,000 

h, m (ft) 

0.11 (0.01) 

0.0011 (0.0001) 



100 

1000 

^x’3 

1000 

10,000 


The results of the estimation process are presented as nondimensional deriva- 
tives . Also presented are the uncertainty levels for each data point . The 
uncertainty levels are proportional to the Cramer-Rao bound and provide an indica- 
tion of the relative uncertainty of the derivatives . Additional information on the 
uncertainty levels can be found in reference 5 . The velocity derivatives were non- 
dimensionalized with respect to the nondimensional velocity V/a (that is , Mach 
number) , and the altitude derivatives were nondimensionalized with respect to the 
nondimensional altitude h/h^. All the nondimensionalizing parameters (q, Vq, a^, 


hp. W. 


and Jy) were average values for each time history . 


In nondimensionalizing the phugoid derivatives , no attempt was made to 
separate the coefficient changes with respect to velocity and altitude from the effects 
due to dynamic pressure . Thus , the nondimensional Z-force derivatives , for 
instance, were related to the dimensional derivatives as follows: 


Z 


V 




z 


h 




instead of by the following, more complete expressions: 






aomVo 



2^S_c 

mVo^ ^0 
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The short-period maneuvers were analyzed at a sampling rate of 10 points per 
second . A study of the effect of sampling rate on the phugoid derivatives is 
summarized in appendix A. The study included sampling rates from 0.1 to 10 sam- 
ples per second , and a sampling rate of 2 points per second was selected on the 
basis of derivative uncertainty levels and computer computation time. Since the 
data were available at 10 samples per second , the subsequent analyses were made 
with 5-point averages of the original data to obtain the 2-sample-per-second data. 


RESULTS AND DISCUSSION 


Short-Period Characteristics 

Eight maneuvers were analyzed to determine the airplane's short-period char- 
acteristics . The derivatives were obtained from elevator pulse maneuvers with 
average bypass door position used as an additional control parameter to represent 
the inlets for the automatic inlet cases . The nondimensional coefficients are shown 
in figure 5 and tabulated in table 1 , with the exception of the inlet derivatives , 
which are shown in a later section . The uncertainty level for each data point is 
also shown . An average value from all eight maneuvers was found for each deriva- 
tive for use in determining the phugoid derivatives . The technique used to find 
the average , which is described in appendix B , weighted each data point on the 
basis of its uncertainty level . Although it was possible to identify the inlet deriva- 
tives in the short-period maneuvers, no significant variations were observed in the 
other short-period derivatives if the inlet model was not included . 

It was not possible to identify the derivatives Z^, X^, , and Xg , and they 

were kept fixed at zero during the analysis . A typical comparison between a flight 
time history and a computed time history is shown in figure 6 . The computed time 
history used the average values of the derivatives , the flight-measured control 
input, and the variable bias terms. The figure shows that the average values 
provide a reasonable representation of the aircraft's flight characteristics. 


Phugoid Characteristics 

Three phugoid maneuvers were analyzed to obtain velocity and altitude deriva- 
tives for the basic aircraft. Cases A and B included the inlet as a separate 
controller . The inlets were fixed for the case C maneuver , so the basic airplane 
derivatives were obtained without having to estimate the inlet derivatives. As 
shown in figure 7 , the fit of the computed and measured time histories is reasonably 
good . The velocity and altitude derivatives and the weighted averages based on the 
uncertainty levels are given in figure 8 and table 2 (a) for the three cases . There 
is a considerable amount of scatter in the data, but, in general, cases A and B agree 
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with case C , validating the use of bypass door deflection as representative of the 
total inlet system . 

Maneuvers A and B were also analyzed to determine the effective velocity and 
altitude derivatives for the combined airplane/inlet system . The computed time 
histories showed essentially the same fit as shown in figure 7 for the basic airplane . 
The resulting derivatives are summarized in figure 8 and table 2(b) and are 
denoted as cases A and B. The differences between the derivatives for the basic 
airplane and the airplane/inlet system are due to the effects of the inlets. Figure 8 
shows that the effects of the inlets on the Z-axis derivatives (C'„ and C' .) as well 

as on the pitching moment derivatives (C' and C' are small. As would be 

expected , the primary effect of the inlets is along the X-axis (which is approxi- 
mately the negative thrust axis) . The average value of C' for the basic airplane 

(fig . 8) is negative , indicating a stalnlizing_tendency with respect to speed changes . 
For the airplane/inlet system cases (A and B, fig. 8) , C' is positive, which would 
indicate a tendency for speed to diverge . '^M 

The uncertainty levels indicate that , in general , cases C , A, and B, in which 
the inlet derivatives were not identified , produced better estimates of the deriva- 
tives . Cases A and B were poorly conditioned for purposes of inlet derivative 
determination , since the bypass doors moved approximately with velocity , causing 
a nearly linear dependence between the velocity and bypass door derivatives. Some 
independent motion did occur because of small atmospheric disturbances , and this 
allowed t^ estimates to be made; but the uncertainty levels are greater than in 
cases C , A, and B. 

In addition to the uncertainty in the phugoid derivatives that resulted from the 
conditioning of the phugoid maneuvers , some uncertainty in the phugoid derivatives 
resulted from the fact that short-period derivatives were kept fixed during the 
analysis of the phugoid maneuvers . The effect of each short-period derivative was 
evaluated by rematching one maneuver (case A) with the short-period derivative 
increased by the value of its average uncertainty level (shown in table 1) . The 
resulting uncertainty levels_are compared with those due to the original phugoid 
time history match for case A in table 3 . The primary short-period contribution to 
the uncertainty in C' and C'. was from C„ . For CL and CL . the primary 

contribution was from C^ and for C and C . it was from C . The table shows 

that the total short-period contribution to the uncertainty of the phugoid derivatives 
is small compared with the uncertainty resulting from the match of the phugoid time 
history, except for the effect of on CL and CL. , where the two contributions 

are about equal . It was assumed that the uncertainties due to the short-period _ 
derivatives in the other cases were the same as those shown in table 3 for case A , 
and when applied to the other cases these uncertainties resulted in the average 
velocity and altitude derivatives and the average uncertainty levels summarized in 
figure 9 . 
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Inlet Derivatives 


Inlet derivatives were obtained from the short-period as well as from the 
phugoid maneuvers, and the results are summarized in figure 10 and table 4. The 
weighted average shown in both the figure and the table is a weighted mean of all 
the points based on the uncertainty levels shown. As a result of the weighting, 
the average values of C „ and C„ were primarily determined from the 

\p \p 

phugoid maneuvers , and the average value of C 
from the short-period maneuvers . 


m. 


was primarily determined 


hp 


During the phugoid maneuvers, the inlet varied primarily with velocity. As a 
means of verifying the consistency of these results with the effects attributed to the 
inlets in the derivatives shown in figure 8 , effective velocity derivatives were 
calculated from the inlet derivatives . A least-squares fit of the bypass door position 
as a function of Mach number yielded a value of equal to -41 percent per 

Mach number. Effective velocity derivatives were calculated using this value. The 
derivatives are compared with the increments from figure 9 in the following table: 



Difference between basic 
airplane derivative and 
airplane/ inlet derivative 
(fig. 9) 

Effective velocity derivatives 
calculated from inlet derivative 
(fig. 10) and 

A6^p/AM^ = -41 percent per M 

AC' 

-0.00259 

-0.0160 




AC' 

0.0526 

0,0462 

AC 

0.00126 

0.00173 

'"m 




Agreement is good for AC' , which is the derivative most affected by the bypass 
doors. The values of AC' are in reasonably good agreement. The agreement of 
AC' is relatively poor, but the magnitudes are low compared with C' , and the 
differences are not significant . 


Aircraft Stability Characteristics 

Eigenvalues were calculated with the derivatives of each case and with the 
average sets of derivatives . All the derivatives were dimensionalized to the same 
flight condition (M = 2.9, h = 21,900 m (72,000 ft)) . One set of short-period 
derivatives, which produced short-period roots of -0.17 ± 1.16/, was used. The 
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phugoid and height mode roots are shown in figure 11. Agreement is reasonably 
good between the roots for the individual cases and the roots from the average 
derivatives , with the exception of case B , the derivatives for which had large 
uncertainty levels (fig . 8) . For the average basic airplane data , the phugoid is 
neutrally damped, with a period of 151 seconds. The height mode is stable, with 
a time to half amplitude of 99 seconds . For the average airplane/inlet system data , 
the phugoid is slightly divergent , with a time to double amplitude of 490 seconds 
and a period of 137 seconds . The height mode is divergent, with a time to double 
amplitude of 114 seconds. 

The airplane is normally flown with a longitudinal stability augmentation system 
in operation. The characteristics of the phugoid and height mode roots for a linear 
representation of the augmentation system are shown in figure 12 . Also presented 
are the specifications for the phugoid mode from reference 6 . It is apparent that 
the augmentation system has little effect on the height mode but improves phugoid 
damping from level 3 to level 1 . However , the low magnitude of the pitch rate in 
the phugoid maneuver with the augmentation system engaged (fig. 7), combined 
with the deadbands within the control system , produced essentially no feedback of 
pitch rate or attitude to the elevator. As a result, the augmented phugoid damping 
is a function of the magnitude of the oscillation . The height mode , not being a 
function of the stability augmentation system , is convergent for the basic airplane 
and divergent for the airplane/inlet system . Limited pilot evaluations of the 
augmented airplane with the automatic inlets (ref. 1) resulted in satisfactory 
(level 1) pilot ratings for the altitude hold task under stable atmospheric conditions 
but unsatisfactory (level 2) pilot ratings when atmospheric disturbances were 
present . 

The effects of the velocity and altitude derivatives on the phugoid and height 
mode roots are shown in figure 13 . The derivatives were varied about the values 
for the average airplane/inlet system . Also shown are the values for the average 
basic airplane . The figure gives an indication of the significance of the differences 
between the derivatives for the fixed and automatic inlets . The Z-axis derivatives 
(fig. 13(a)) show little effect of the inlets on the phugoid and height modes. The 
X-axis derivatives (fig. 13(b)) show large differences between the basic airplane 
and the airplane/inlet system derivatives , but X^ has little effect on either the 

phugoid or the height mode . The derivative Xy (a positive value indicating an 

increase in net propulsive force with speed) is shown to be a primary derivative for 
the height mode root. The pitching-moment derivatives are shown in figure 13(c) . 
The differences in My due to the inlet configuration are large but have little effect 

on the roots . The differences in have a primary effect on the phugoid frequency . 

Thus , when the basic airplane and the airplane/inlet system derivatives in figure 9 
are compared, only the differences in CL and C' ^ appear to be significant. 
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CONCLUDING REMARKS 


The longitudinal stability and control derivatives of the YF-12 airplane were 
determined from flight tests at a Mach number of approximately 2 . 9 by using a 
maximum likelihood parameter estimation technique . Short-period pulses and 
phugoid maneuvers were analyzed to obtain the velocity and altitude derivatives 
as well as the short-period derivatives . Engine inlet bypass door position was 
successfully used to represent the total inlet system , and the effect of the inlets on 
the velocity and altitude derivatives was determined . 

The phugoid mode of the basic airplane had neutral damping , and the height 
mode was stable . With the addition of inlets in the automatic configuration , the 
phugoid mode was slightly divergent and the height mode was divergent with a 
time to double amplitude of approximately 114 seconds. The results of the deriva- 
tive estimation indicated that the change in the height mode characteristics was 
primarily the result of the change in the longitudinal force derivative with respect 
to velocity . The net propulsive force of the basic airplane decreased with increasing 
velocity , whereas the net propulsive force of the airplane/inlet system increased 
with increasing velocity . 


Dryden Flight Research Center 

National Aeronautics and Space Administration 
Edwards, Calif. , January 31, 1977 
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APPENDIX A. -REQUIREMENTS FOR TIME HISTORY 
SAMPLING RATE AND LENGTH 


In the determination of stability and control derivatives , there is a tradeoff 
between the sampling rate (or computation time) and the accuracy of the results . 
The tradeoff between computation time and accuracy was evaluated for the phugoid 
maneuver denoted case A . Data acquired at a sampling rate of 10 points per second 
were available. Sampling rates down to 0.1 point per second were obtained by 
thinning the original time history . The identification algorithm was iterated until 
the cost functional did not change more than 0 . 1 percent . 

The velocity and altitude derivatives and their uncertainty levels are shown in 
figure 14 as a function of sampling rate . Agreement between the 10-point-per- 
second data and the lower sampling rate data is good for sampling rates down to 
1.25 points per second. Below this rate, agreement becomes poor, and at a 
sampling rate of 0 . 1 point per second the estimation method does not converge on a 
solution. The uncertainty level increases as sampling rate decreases; however, 
this may be due to the noise characteristics of the data rather than a true indication 
of accuracy . 

For data acquired at sampling rates of 1.25 points per second and above, con- 
vergence was obtained in six iterations. At sampling rates of 0.625 and 0.25 point 
per second, seven iterations were required for convergence. A sampling rate of 
2 points per second was chosen as a reasonable compromise between computation 
time and accuracy . Before the final analysis , the 10-point-per-second data were 
presmoothed using a simple 5-point average to obtain the 2-point-per-second data. 

Another factor that affects computation time and the accuracy of the results is 
the length of the time history (or observation time) from which the derivatives are 
estimated. The phugoid maneuver of case A (with 2-point-per-second data) was 
used to evaluate the effects of observation time on the velocity and altitude deriva- 
tives, and the results are shown in figure 15. Since there were no control inputs 
and the oscillations were nearly uniform , all the time histories had the same 
starting point . The derivatives shown for observation time T correspond to the 
altitude time history shown in the figure from 0 to T . Although estimates of the 
derivatives were obtained for as little as one-half cycle of the phugoid oscillation 
(T = 60 seconds) , approximately two cycles of oscillation appear to be necessary 
before reasonable estimates of all the derivatives can be obtained. 
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APPENDIX B . -USE OF UNCERTAINTY LEVELS IN THE CALCULATION 


OF AVERAGE DERIVATIVE VALUES 


When several estimates of a derivative are available , it is desirable to use 
information about the relative accuracy of the estimates to obtain the best overall 
estimate of the derivative . The parameter estimation method in reference 4 provides 
an indication known as uncertainty level of the accuracy with which the derivatives 
have been identified . Uncertainty level is analogous to standard deviation , and 
further discussion of uncert^nty level can be found in reference 5 . For this study , 
an average derivative value d was determined by using a weighted average as 
follows: 


N 


d = ^-=1 


W.d. 

i i 


N 


2 


i=l 


where 


W. 


N 


ith estimate of the derivative 
weighting for the ith estimate 
number of estimates 


The expression for the weighting was based on uncertainty level and was given by 

_2 

W. = ^Uncertainty level of d.^ 

An average uncertainty level u was determined from the weighted average 

N 


2 *i“(‘ 


-2 i=l 
u = 


N 


2 


i=l 


where is the uncertainty of derivative d^. For the weightings used in this study, 


APPENDIX B. -Concluded 


this expression reduced to 
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TABLE 1. -SHORT-PERIOD DERIVATIVES 


Case 

Value of — 

per deg 

per deg 

per deg 

S’ 

per rad 

* 

TW ^ 

per deg 

A1 

Short-period 

derivative 

Uncertainty 

level 

-0.028960 

±0,004770 

0.0008532 

±0.0021500 

-0.0009680 

±0.0000109 

-1.4510 

±0.2511 

-0.0013100 

±0.0000658 

A2 

Short-period 

derivative 

-0.031510 

0.0012890 

-0.0009712 

-0.9344 

-0.0011500 

Uncertainty 

level 

±0.002868 

±0.0018405 

±0.0000133 

±0,1555 

±0.0000378 

B1 

Short-period 

derivative 

-0.026610 

0.0015230 

-0.0008586 

-0.9145 

-0.0008808 

Uncertainty 

level 

±0.004756 

±0.0019235 

±0.0000093 

±0.1737 

±0.0000392 

B2 

Short-period 

derivative 

-0.037660 

0.0018830 

-0.0008917 

-0.6898 

-0.0009513 

Uncertainty 

level 

±0,011280 

±0.0044365 

±0.0000151 

±0.3874 

±0.0000718 

Cl 

Short-period 

derivative 

-0.026640 

0.0004785 

-0.0008647 

-1.4520 

-0.0010350 

Uncertainty 

level 

±0.002439 

±0.0019815 

±0.0000171 

±0.2494 

±0.0000485 

C2 

Short-period 

derivative 

-0.031470 

-0.0000476 

-0.0008802 

-1.0530 

-0.0009537 

Uncertainty 

level 

±0.003916 

±0.0028345 

±0.0000207 

±0.2579 

±0.0000484 

T\ 

Short-period 

derivative 

-0.043340 

0.0013520 

-0.0009257 

-0.9019 

-0.0011900 

U 

Uncertainty 

level 

±0.004632 

±0.0020700 

±0.0000123 

±0.2174 

±0.0000525 

g 

Short-period 

derivative 

-0.022660 

0.0008864 

-0.0009715 

-1.3600 

-0.0010180 


Uncertainty 

level 

±0.004797 

±0.0020670 

±0.0000116 

±0.1951 

±0.0000422 

Weighted 

Short-period 

derivative 

-0.029782 

0.0010191 

-0.0009200 

-1,0853 

-0.0010462 

average 

Uncertainty 

level 

±0.003845 

±0.0021827 

±0,0000127 

±0.2134 

±0.0000474 
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TABLE 2.-PHUG0ID DERIVATIVES 
(a) Basic airplane 


Case 

Value of — 

r' 

^ 7 


C 

a A 

CL 


A 

Phugoid 

derivative 

-0.0590280 

0.3209103 

-0.0341026 

0.0358694 

-0.0002269 

-0.0069156 


Uncertainty 

level 

±0.0202356 

±0.0432446 

±0.0069617 

±0.0122550 

±0.0006039 

±0.0013324 

B 

Phugoid 

derivative 

-0.0051690 

0.3621604 

0.0162794 

-0.0529985 

0.0042093 

-0.0025460 

Uncertainty 

level 

±0.0352015 

±0,0460856 

±0.0141296 

±0.0143139 

±0.0014389 

±0.0018926 

n 

Phugoid 

derivative 

-0.0682926 

0.2718828 

-0.0276689 

-0.0072000 

0.0021633 

-0.0008271 


Uncertainty 

level 

±0.0359928 

±0.0331589 

±0.0028062 

±0.0041600 

±0.0010258 

±0.0009657 

Weighted 

Phugoid 

derivative 

-0.0499972 

0.3077680 

-0.0270963 

-0.0062877 

0.0008304 

-0.0028689 

average 

Uncertainty 

level 

±0.0273144 

±0.0395794 

±0.0044335 

±0.0065784 

±0.0008477 

±0.0012517 


(b) Airplane/inlet system 


Case 

Value of — 

^z 

r' 

r' 

X 

M 

CL A 

CL 

’^M 

C' A 

m* 


Phugoid 

derivative 

-0.0498299 

0.3185730 

0.0177792 

-0.0716930 

-0.0008564 

-0.0060487 


Uncertainty 

level 

±0.0106757 

±0.0375977 

±0.0017023 

±0.0088026 

±0.0003428 

±0.0011962 

B 

Phugoid 

derivative 

-0.0731338 

0.4088526 

0.0417087 

-0.0406737 

0.0048611 

-0.0037563 

Uncertainty 

level 

±0.0301993 

±0.0504131 

±0.0024658 

±0.0066243 

±0.0012091 

±0.0020923 

Weighted 

Phugoid 

derivative 

-0.0524187 

0.3508400 

0.0255026 

-0.0518891 

-0.0004309 

-0.0054840 

average 

Uncertainty 

level 

±0.0142344 

±0.0426229 

±0.0019811 

±0.0074854 

±0.0004664 

±0.0014686 



























TABLE 3. -EFFECT OF SHORT-PERIOD DERIVATIVE UNCERTAINTIES ON PHUGOID DERIVATIVES 



Uncertainty in 
short-period 
derivative , 



Phugoid derivative 



Short-period 

derivative 


c* 

CL 



m* 


percent 


Uncertainty in phugoid derivative due to 
short-period derivative uncertainty, percent 



12.9 

4.42 

2.76 

0.36 

7.00 

2.69 

0.17 

a 

214.1 

0,15 

8.02 

0.33 

0.13 

16.92 

0.15 

c 

"’a 

1.4 

0.08 

0.31 

1.00 

0 

0.29 

1.22 

c 

19.7 

1.36 

0.14 

2.81 

0,66 

0.44 

2.23 

e 

4.5 

0.17 

0.81 

1.92 

0.03 

0.64 

0.03 

Total uncertainty due to 
short-period derivatives , 







^ ^ (Uncertainty ) 

4.63 

8,53 

3.58 

7.03 

17.16 

2.55 

Uncertainty due to 
phugoid time history match 
for case A 

21.42 

9.57 

40.03 

11.80 

12.28 

19.78 

Total uncertainty 
for case A 

21.91 

12.82 

40.19 

13.74 

21.10 

19.94 
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TABLE 4. -INLET DERIVATIVES 




Inlet derivatives and 




uncertainty levels 

Case 

Value of — 

bp 

C y t 

^bp 




per percent 

per percent 

per percent 

A 

Derivative 

Uncertainty 

-0.0000457 

-0.0021895 

0.0000321 


level 

±0.0005731 

±0.0002774 

±0.0000215 

B 

Derivative 

Uncertainty 

0.0012055 

-0.0006501 

0.0000019 


level 

±0.0005095 

±0.0002241 

±0.0000196 

A1 

Derivative 

Uncertainty 

0.0003262 

-0.0004757 

0.0000998 


level 

±0.0021275 

±0.0012145 

±0.0000168 

A2 

Derivative 

Uncertainty 

0.0003846 

-0.0005902 

0.0000467 


level 

±0.0013875 

±0.0011160 

±0.0000084 

B1 

Derivative 

Uncertainty 

-0.0003080 

-0.0003555 

0.0000387 


level 

±0.0006220 

±0.0004877 

±0.0000055 

B2 

Derivative 

Uncertainty 

0.0002862 

-0.0001223 

0.0000203 


level 

±0.0038900 

±0.0016910 

±0.0000302 

Weighted 

Derivative 

Uncertainty 

0.0003887 

-0.0011230 

0.0000422 

average 

level 

±0.0007637 

±0.0003925 

±0.0000102 


E-24984 


Figure 1. YF -12 airplane. 









Figure 3. Schematic of inlet control system. 



Figure 4. Axis system. Arrows indicate positive directions of values. 
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(a) Phugoid case A . 

Figure 7. Comparison of flight and computed time histories 
for phugoid maneuvers . 


®bp’ 

percent 


'e« 

deg 


X’ 

9 


n’ 

g 


10 p 

0 


Flight 

Computed 


-10 
-5 
-6 
-7 

. 20 ,- 

.10 


J I I I 


J 1 I I ^ I 


0 

1.40 

1.00 

.60 


_l I ^ I L_ 


_l ^ I ^ L_1 



-| 77.000 

- 74,000 

- 71,000 


h, 

ft 


9 , 

deg/sec 


V, 

m/sec 


a, 

deg 


.40, 


-.40 




1 ^ ^ L 



3000 
2900 
-• 2800 


V, 

ft/sec 


0 40 80 120 160 200 240 

t, sec 

(b) Phugoid case B. 
Figure 7. Continued. 
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Figure 8. Velocity and altitude derivatives obtained from 
phugoid maneuvers , 
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Figure 9. Summary of velocity and altitude derivatives and average 
uncertainties for M 2.9, h 21,900 m (72, 000 ft) . 




Figure 11. Phugoid and height mode roots for individual cases and weighted 
average derivatives from figure 8 for M = 2.9 and h = 21,900 m (72, 000 ft) . 


35 







Derivative 

Z„ 


Basic 

airplane 

value 

o 

A 


(a) Zy and Z^. 

Figure 13. Phugoid and height mode characteristics as a function of velocity and 
altitude derivative values . 
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Figure 13. Concluded 
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Figure 14. Effect of_sampli 
phugoid maneuver A . 
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Figure 15. Effect of observation time on velocity and altitude derivatives for 
fixed sample rate of 2 points per second with phugoid maneuver A . (Derivatives 
shown at T correspond to observation of time history from 0 to T .) 
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